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1. INTRODUCTION 


We consider discrete approximations to the scalar conservation law 

(1*1) |-jr [u(x,t)] [f (u(x, t) ) ] = 0, (x,t)€Rx[0,»). 

Let v(t) = {v^(t)} be the approximate solution, and denote by 

(1.2) TV[v(t) ] = l |Av v+1/2 |, Av v+1/2 = v v+1 (t) - v v (t), 

its total-variation at time level t. A desirable property for such an 
approximate solution to share with the exact one, is that its total variation 
should decrease in time*. Difference scheme which give rise to such total- 
variation diminishing solutions — called TVD schemes after Harten [3] — is the 
subject of this paper. 

TVD schemes prevent spurious oscillations in their solutions, and unlike 
monotone schemes, they can still allow for high-accuracy in most of the 
computational domain. Consequently, the TVD schemes can offer a substantial 
gain in computational efficiency as indeed was verified in a wide range of 
applications, e.g. [11] and the references therein. 

Sufficient TVD criteria for explicit and implicit fully-discrete schemes 
were given by Harten in [3], [4], and analogously for semi-discrete schemes in 
[13], [7]. Necessity for three-point schemes was proved in [15, Lemma 2.2] 
and a general TVD characterization for multi-point stencils was provided in 


^We use the notion of order in its weak sense; thus, decrease means 
nonincrease, positive refers to nonnegative, etc. 
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[5], [9]. Roughly speaking, these criteria assert that a given scheme has the 
TVD property, provided it can be written in an appropriate incremental form 
which meets a certain positivity condition, augmented with a CFL restriction 
in the explicit case. A difference approximation of (1.1) can be equally 
represented by a variety of different incremental forms, yet the positivity 
condition mentioned above is not invariant under such different 
representations. Hence, the key step in seeking the TVD property according to 
the above criteria, requires us to find the correct incremental form so that 
this positivity condition could be met, e.g. [10, Part I]. 

In this paper, we provide alternative more convenient TVD 
characterizations, in the sense that they are uniformly valid for the various 
different incremental representations of a given scheme. To this end, we 
first note that the total variation of a grid function depends solely on its 
extrema values, see (1.2). It is therefore plausible to assert that in order 
for a difference scheme to share the TVD property, its incremental 
coefficients should be controlled only at critical neighborhoods where the 
approximate grid solution attains extrema values. Indeed, our sufficient TVD 
conditions have the flavor of this assertion, namely, they place a positivity 
restriction only on those incremental coefficients which are associated with 
such critical neighborhoods. Moreover, we are also able to express these 
local TVD conditions solely in terms of the numerical fluxes of the nonlinear 
schemes, rather than invoking any of their special incremental 
decompositions. Putting it in different words, we show that the TVD property 
holds for difference schemes whose numerical viscosity correspond to upwind 
differencing at extrema values but can be arbitrary otherwise. Thus, in 
contrast to the more restrictive global positivity conditions mentioned 
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earlier, our TVD criteria are localized to extrema values, and consequently 
can be equally applied to different incremental representations. 

We begin by discussing the semi-discrete case in Section 2. Fully- 
discrete implicit and explicit schemes are treated in Section 3. To utilize 
our TVD criteria in the latter cases, the extrema values at the next time- 
level are to be known in advance. This necessitates additional ingredients in 
the fully-discrete cases, whose purpose is to provide us with such a priori 
control on the behavior of extrema at the next time level. In particular, 
standard recipes of constructing TVD schemes which make use of anti-diffusive 
correctors and limiters, are all shown to naturally follow in light of our 
above arguments. 


2. SEMI-DISCRETE SCHEMES 

We consider semi-discrete schemes in the conservative form 

(2.1) ~ — v (t) = — tt— — [h i< — h 

dt v Ax y v+ 72 v- 72 

with Ax = 4 (x , - x ,) being the variable meshsize and h , 

v 2 v+1 v-1 v+ 72 

denotes the Lipschitz continuous numerical flux which is consistent with the 
differential one 


( 2 . 2 ) 


h i = h(v 
v+ 72 v-p+1 



h(u,u, . . .u) = f (u) . 


To study the TVD property of these schemes, we forward difference (2.1), 
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( 2 . 3 ) Av . i. = - — - — [h ., /0 “ h 0.1/1 + 7 — [h i, - h ,.] , 

dt v+ V2 ^ X v + 1 v + 3/2 v+ V2 Ax y v+V2 v- V2 


multiply by s v+ l^“ s^ + ^(t) = sgn[Av^ + ^(t) ] and sum by parts, obtaining 


(2.4) ^IV[v(t)) -I Is v+1/2 - Vl/2 1 • IVl /2 - VV 2 ] 


The only contributions to the sum on the right came from extrema values 
where s v+ * s v _ 1 / » an< ^ requirement of these contributions to be 

negative yields 


Lemma 2.1: The semi-discrete scheme (2.1) is TVD, if we have 


(2.5a) h 1 . > h 1 . at maxima values v (t), 

v+V 2 — v- a /2 v 


(2.5b) h 1 . < h 1 , at minima values v (t). 

V+V 2 — v- V 2 v 

In other words, Lemma 2.1 requires maxima values to decrease in time and 
minima values to increase in time. Moreover, if the distance between such 
extrema values exceeds the stencil width of 2 p + 1 cells, then the 
corresponding terms inside the summation on the right of (2.4) are independent 
and consequently (2.5) is also necessary for TVD in this case. 

We now turn to discuss the relation between the TVD criteria in Lemma 2.1 
and a different kind of TVD conditions due to Harten [3], [4] and Osher [9]; 
see also [5], [13]. In order to implement the latter, one should start with 
nonlinear semi-discrete schemes which assume the incremental form 
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( 2 . 6 ) 


dt 


, (t) ‘ c v+ 1/2 V v+ 1/2 " V l/2 iV v- 1/2 • 1/2 5 W 0 - v v (t) - 


The nonlinearity is reflected here by the possible dependence of the 
coefficients C* + y on v y _ p+1 , . . . ,v v+p . Forward differencing of (2.6) 
gives 


( 2 . 7 ) |jAv v+1/ - (C 


v+3/2 Av v+3/2 



Multiplying (2.7) by s v+ \j and summing by parts we find 


(2.8) |Av v+1/ j = l s v+l/ - ^ Av v+ ^ = 


-I K 

v 


s v+ V2 S v- ^2 ^v+ V2 + ^ S v+V2 S v+3/2^ %+ I/2 AV v+ V2* 


and using the fact that 2 Av y+ ^ = s y+ ^ | A Vy+ 1; J where S v+ ^ = 1, we end 
up with 


(2.9) ^TV[v(t)] - - HU - Vl/ 2 Vl/2 ) ' C v + l/2 + 

(1 ' Vl/ 2 % + 3/2 >C v + l/ 2 1-|4v v+l/2 | - 

The quantities inside the two parenthesis on the right are positive. Hence, 
the summation on the right is positive and consequently the scheme (2.6) is 


2 


The signum function at zero is defined to be 


±1 , so that its square = 1, 
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TVD, provided the incremental coefficients, C* + \j , are positive 

< 2 - 10 > < + v 2 >°. 

The positivity requirement (2.10) is the usual condition which characterizes 
the TVD schemes (2.6), e.g., [3], [5], [7], [13]. Given a semi-discrete 

conservative approximation of (1.1), it can be equally represented in a 
variety of different incremental forms. The positivity condition is not 
invariant, however, under such different representations. Thus, a key step in 
seeking the TVD property for a given scheme, requires us to find the correct 
incremental form so that the positivity condition (2.10) could be met, e.g. 
[10], 

Lemma 2.1 provides us with a local TVD criterion which makes no reference 

to the incremental representation of the scheme (2.1). How does this compare 

with the global positivity condition placed on the incremental coefficients in 

(2.10)? a second glance at (2.9) shows that whenever the grid values v^ 

and v v+ i are located in a monotone profile, i.e., when both 

1 - s v _ 1^2 s v+ I /2 an< * 1 - s v+ l/ 2 s v +3/2 vanish, then the corresponding term 

in the summation on the right of (2.9) also vanishes independently of the 

incremental coefficients C* 1 . . This tells us, therefore, that the 

v+ V 2 

positivity condition (2.10) can be localized to extrema values, bearing a 
close similarity to the local nature of the TVD criterion (2.5). To be more 
precise, let us abbreviate 
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then (2.9) reads 


(2.12) ^ TV[ v(t) ) - - l U V C l. l/2 + X v+1 C; +1/2 ] • |AV v+1/2 | . 

and we are led to the following. 


Lemma 2.2: The semi -discrete scheme (2.6) is TVD, if we have 


(2.13) 


x v^v+l/2 + ^v+l^v+l/ 2 — ■ 


For smooth grid functions, we have almost everywhere (i.e., with the exception 
of critical neighborhoods), x v = X v+ j = 0» hence the TVD condition (2.13) is 
automatically fulfilled in these cases. 

Once the positivity condition (2.10) was localized to those incremental 
coefficients associated with extrema values (2.13), we can go one step further 
and complete the comparison with Lemma 2.1, dealing with the numerical fluxes 
instead. 

To this end, the scheme (2.1) is rewritten in its canonical incremental 
representation (2.6) where, see [15, Section 2] 


(2.14) C 


f ( v ) - h 


v+ V 2 Ax 


A v 


v+ l h 


v+ 1/2 


’ ^v- 1/2 Ax 


f ( v ) - h 


4 V V 2 


v- V 2 


, Av 1/ ^ 0. 
v+ 72 


Applying Lemma 2.2 to these coefficients, then (2.13) reads 


»V f(T V 5 ‘ h V+ Vi. X V+1 f( '»H ) ' VV? , n 

hr • 7, SiTTT ' 4V — z ^ 


V+ I/2 


v+1 


v+ l/ 2 


(2.15a) 
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or, equivalently, 


«-‘5M s v+1/ (h v+1 - f(v,» +T ^i- . (h_ - f(v_, » <0. 


v '2 v+ '2 - ““v+1 


IT— • V VW 2 _ IV Vi ;; i 


Hence, in case of an isolated extrema value where x y = 2 , x v±1 = °* the 
inequality (2.15) is fulfilled if and only if 


(2.16a) 


V V 2 <h v+ 1/2 " - °’ 


(2. 16b> 


Vi/ 2 (h »-v 2 ' £<v '> )) - °- 


In fact, (2.16) covers the general case of extrema values whether 

isolated or not. For, if x = x , , = 2, then, since v is also an 

’ A v v+1 v+1 

extrema value, we have in view of (2.16b) 


(2.17) 


s i .(h i f(v )) <0, 
v+ V 2 v+ V 2 v+1 — 


and a weighted average of (2.16a), (2.17) yields (2.15). 

We summarize this by stating 

Corollary 2.3: The semi-discrete scheme (2.1) is TVD, if we have 


(2.18a) 


h I, > f(v ) > h l, at maxima values 

v+ a / 2— v — v- V 2 


v v ( t), 


(2.18b) h i. < f(v ) < h l, at minima values 

v+7 2 — v — v- V 2 


v v (t). 
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The TVD condition (2.18), which was derived on the basis of the 
incremental decomposition (2.14), is somewhat more stringent than our TVD 
criterion (2.5) in that the former requires f(v^) to separate between the 
numerical fluxes on both sides of extrema values. Incremental decompositions 
of (2.1) other than (2.14), may lead to slightly different local TVD 
conditions; yet, they all share a similar kind of a separation requirement at 
extrema values, which in view of the consistency relation (2.2) is a generic 
property of the TVD numerical fluxes. 

Lemma 2.1 and Corollary 2.3 enable one to verify the TVD property of 
first as well as higher order accurate semi-discrete schemes, without making 
reference to any of their special incremental representations. To demonstrate 
this point, we turn to 


Example 2.4: Consider the class of generalized MUSCL schemes [8], where 


(2.19) 


, ,E, . Ax , Ax , 

h v+v 2 " h <v v + rV Vh - — Vi>' 


Here, Ax^ = Ax is the uniform mesh spacing, h (•,•) stands for any 
E-flux, satisfying 


( 2 . 20 ) 


sgn(w y + 1 - w y ) • ( h ( w v » W v+1^ " f(w)) _< 0 


for all w between w^ and w v+ j> and d y is an approximate derivative at 
x^ which guarantee second-order accuracy if chosen so that 


( 2 . 21 ) 


Ad i. = -i— • [Av i. - Av ij + 0 (Ax)" 
v± V2 Ax v+ 1 2 v ~ ' 2 
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In [8, Lemma 2.3], Osher introduces a special incremental decomposition 
of these schemes in order to show that they meet the positivity condition 
(2.10) and hence share the TVD property, provided for each v we have 


(2.22) 0 < . ■ • d <1. 

— Av , i. v — 
v± V2 

Note that in the particular case of v^ being an extrema value, (2.22) 
implies that d^ must vanish and, consequently, that accuracy degenerates to 
first-order at these points. 

In contrast to the special positivity arguments made above. Lemma 2.1 
suggests a straightforward TVD derivation in this case. Localized at extrema 
values, we set d^ = 0, so that in view of the E-condition (2.20), TVD is 
guaranteed if in addition we have 


(2.23a) 


sgn(v 


v+1 


V = s S n(v v+l 


- v 


Ax , v 

tW 


(2.23b) 


sgn(v v - v yl ) = sgn(v 


Ax , % 

v " V v-1 2 d v-l )s 


i.e., if the neighboring discrete derivatives of extrema values satisfy 


(2.24) 


1 d v+i 1 1 K 

v+ l 2 


V2 


One possible choice for such discrete derivatives, d^, could be 


(2 * 25) d v = S7- Min MAv v _i /2 |, |Av v+1/ J], s v = i[s v _ 1/2 + s v+1/ ^ 
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The above TVD analysis relies on the conservative form of nonlinear 
difference schemes. We turn now to discuss another representation which is 
found useful to utilize TVD criteria for such schemes, making use of their 
viscosity form. To this end, recall the definition of the incremental 
coefficients, C* + jy, in (2.14). The identity 


(2.26) Ax ..C . i, - Ax C + . = - V+ l/ l 

v+1 v+ V 2 v v+ V2 Av^ + i j 


Af 


v+ V 2 


= " f(v v ) ’ 


shows that between these two incremental coefficients, C* + ^ , there is only 
one degree of freedom, which could be expressed in terms of Q^ + , 


(2.27) 


Q , 1 . = Ax ,C 1 . + Ax C 1 .. 
v+ V2 v+1 v+ V2 v v+V2 


Eliminating C, 1 . from (2.26) and (2.27) we find 
v+ V 2 


(2 * 28) Cv+ V2 2Ax v+1 (Q v+ V2 + Av v+ 1/2 J ’ C v+ V2 2Ax v (Q v+ V2 Av^ + ^ ) ' 


Af 


Our scheme (2.1) is then recast into the viscosity form [16] 


(2.29) ^v v (t)- 




«Vi>l + 


+ 24 ^ ‘V !/£ V V 2 ' V 1 / 2 4 V ‘/ 2 1 1 


thus revealing the role Q plays as the numerical viscosity coefficient, 
Applying (2.13) to C* + y § iven in ( 2 * 28 )« then Lemma 2.2 tells us 
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Lenma 2.5: The scheme (2.29) is TVD, if its viscosity satisfies 


Af 


(2.30) 


(x Ax . + x , ,Ax ) Q 
V v+1 A v+1 v 


V+V 2 — (X v AX v+l X v+l““v' Av 


, Ax ) 


v+ V 2 


v+ V 2 


In the case of equally spaced meshpoints, Ax^ = Ax, we conclude that the 
scheme 

<2 - 3I) If v v (t> ' ‘ ih tfCv v+i 5 ‘ «V|H + 

+ JKZ [Q v+ I /2 V v+ V 2 " V V 2 Av v- V2 1 ’ 

is TVD, provided the following simple Inequality is fulfilled at neighborhood 
of extrema values 

< 2 ' 32) ( *v + WVl/ 2 i (x v ‘ W • 

Example 2.6: Consider a first-order accurate TVD scheme of the form 


<2 - 33) If v v (t> ■ - 4^ [ V l/ 2 " V 2 ' ' 

■ z hr- 1£<v v.i ) ' «vi>i + 2 sr [ Vv2 v v+V 2 ' Vu'Vn 1 ' 

In order to convert it into a second-order accurate TVD scheme, we add to it 


an anti-diffusive conservative difference 


-13- 


(2.34a) ~ v (t) = - [h . - h ,.] - -±- [h, - h ,.]. 

dt v Ax 1 v+V2 v-V2 Ax v+V2 v- 72 


The numerical flux correction, h i , , is chosen to be of the form [10, 

v+ 72 

Corollary 4.9] 


(2,34b) *\>+l/ 2 2 ^ g v + g v+l V V5&+1 


where the so called modified flux correction, , should satisfy 


(i) sgn(g^) = s^ = s v+ ]y , at nonextrema values ( t ) , 


(ii) 

Ay 

g v 

= o, 



at extrema 

values 

v v (t) 

• 


Now, 

if 

v v (t) 

or W C) 

is 

an extrema 

value , 

then 

by (ii) 

we have 

that 

Ay 

= 0 or 

o 

ii 

i-H 

+ 

j wP 

and 

consequently 

yw 

h . 
v + 

V2 

vanishes 

in both 

cases 

since 

by (i) 









\, + V 2 * 7 (i v + l - S v + 1/J* v+I l) - 0 or h u+ 1/ 2 ■ -j (S v - vvj«vl ] ' °' 


Hence, h i. = h = 0 at extrema values v (t) and in view of Lemma 
v+ V2 v- 72 v 

2.1, the modified scheme (2.34) inherents the TVD property of (2.33). Next we 
observe that the modification of (2.33) into (2.34) has the net effect of 


decreasing the original first-order viscosity, Q^ + , into 

/v ^ 

h_ _ 

accuracy [17, Lemma 4.4], the latter 


h + i/ 

Q , 1.-2 — 2-2- ; for second order 

^v+ 72 Av v+1/2 ’ 

should be a Lipschitz continuous grid function of order 0(|Av^ + ^), i.e.. 
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V 1/2 " I V I /2 V v+ l/ 2 + 0(4V v+ V 2 )2 - 


To this end, one could choose the modified flux correction, , as 


(2.34c) S v - / • Hln[|Q u+1 ^v v+1/2 | > | %_ ^ 5 i [...l/j 4 % + 1/J 


In this way, second-order accuracy is achieved away from extrema values. 


noting that 


>+ V2 


takes there the value 


1 


° r g v + l 


which modulo 


quadratic error terras are equal — Q^ + y^ v+ 1/ , 


Example 2.7 : A simple recipe suggested by (2.32), for constructing a TVD 
scheme with second-order accuracy away from extrema values, is to set the 
numerical viscosity Q to be 


Af + 1/ 

(2 - 35> Vy 2 - ( *v -W 

The resulting scheme (2.31), (2.35), amounts to the usual second-order central 
differencing augmented with first-order conservative connection at extrema 
values 


(2.36) |jv v (t) - - J-[f(» v+1 ) - f( Vl )) - 2n[*X v+ l /2 4f v+ l /2 - 4X V _ V 4£ v _ y 2 ] • 

The last two examples dealt with TVD schemes which are second-order 
accurate away from extrema values. In our final example for this section, we 
demonstrate a simple recipe of enforcing the TVD property on arbitrary 


conservative schemes while maintaining their high accuracy away from extrema 
values, see also [9]. 
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Example 2 . 8 : Let ^ v + 1/ be any highly accurate consistent flux. For 
example, the 2p-th order accurate central differencing are identified with the 
numerical flux, e.g., [18] 

P 

h i. = I d, [f(v ,) + ... + f(v . )] , 
v+ V2 k ti k P v+k v-k+1 ’ 


where 


a- 1 , , _ 2 a- 1 

d n 2 and d 12 3 » d 22 12 


for second and fourth order accuracy, or 


A _ (-1) Ax , _ 2n 

kN _ . .kAxv ’ = 2N+1 

2sin(— ) 


for spectral accuracy occupying periodic stencils of 2N+1 meshpoints. 
Next, we denote by 


(2.37) 


h v+V 2 = I ! «V + f(, vtl> - W 2 ‘ |f<, v+l> ' £<v v ) l ] 


the usual first-order accurate upwind numerical flux, and let us consider the 
semi-discrete conservative scheme 


(2.38a) 


d 

v 

dt v 


(t) = - J~ [H i, - H i.], 
Ax y v+ V2 v- /2 


where H 


v+ V2 


is defined as follows 
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(2.38b) 


^v+ V2 |s vVi ! *\>+V2 + ^ | 8 vVl |)h v+V 


s v 5 2 <V V 2 + V *A?" 


Away from extrema values |s^| = |s^ )+1 | - 1> and the original high accuracy 

of V v 2 - \ + v 2 


V ' ' Vdel 1 

is retained in those regions. At extrema values 


3 = 0, hence H , 1 , coincides with the upwind flux h U , satisfying 

v v+ V2 v+ '2 


h" + l /2 - h v- l/ 2 " 2Ma *l f (v v ), «v v+1 )] - 2Min|f(v v . 1 ), f(v v >] > 0, 


at maxima values, and the inverse inequality 


h l + l/ 2 ~ h v- l/ 2 = 2Mln t f ( v v ). f < v v+ i>3 " 2Max[f(v v _ 1 ), f(v v >] < 0, 

at minima values. Consequently, the scheme (2.38) is TVD by Lemma 2.1. 

We note that the numerical flux in (2.38b) is in general not 

smooth, except for the second-order case, p = 1, where the scheme (2.38) 
coincides with the previous example (2.36) and its global second order 
accuracy is maintained, e.g., [3], [4], [9]. The highly accurate stencils, 
p > 1, require further numerical and analytical investigation with regard to 
their accumulated accuracy in extrema free regions. 


Remarks: (i) It is instructive to see why the necessary and sufficient 
conservative TVD criterion in Lemma 2.1 is reduced to the sufficient 
incremental TVD conditions derived from (2.12). To this end, let us insert 
the incremental coefficients (2,14) into (2.12) obtaining 


-17- 


d_ 

dt 


TV[ v(t) ] = 



f(v ) - 
V 


A v 


hi i . 

_v+J/ 2 + 


"v+1 


v+ V 2 


Ax 


v+1 


f ( v , , ) “ 
v+1 


1 1 . 

v+ V2 


Av 


v+ V 2 



Now, Lemma 2.2 and Corollary 2.3 were derived by requiring a termwise 
positivity of the brackets inside the summation on the right, see (2.15). 
Instead, if we first reindex this summation writing it as 


-I 


v 



[s 


V+ V2 


(f(v ) - h i.) + s 1 . ( f ( v ) - 
v v+ V2 v- V 2 v 



we then end up with the necessary and sufficient TVD criterion (2.5). This 
makes apparent the difference between the two derivations due to the 
nonlinearity of the schemes. 

(ii) The inequality (2.32) shows that the scheme (2.31) has the TVD 
property with an arbitrary amount of viscosity, except for intervals 
containing isolated extrema values where we need at least 


(2.39) 



Af 


Av 


y ±}/2 


v+ 72 


The quantity on the right corresponds to upwind differencing, and is 
responsible for the familiar first-order 'clipping' phenomenon at the extrema 
of TVD schemes, e.g., [3], [7], [10]. 


(iii) A classical argument which involves Helly's theorem, Lipschitz 

continuity of lv(»,t)Ll and the diagonal process implies the convergence 
* 1 Li 

of TVD schemes to a weak solution of (1.1), e.g., [4], In particular, this is 
true for central differencing, Q^ + \j = 0, augmented with extrema upwind 


18- 


differencing (2.36). However, the limit solution may still be a physically 
irrelevant one, e.g., [8]. To avoid the latter, say in the convex case 
where f"(u) > 0, it is enough to have viscosity at the amount which exceeds 
[17] 


(2 ' 40) «v + 1/ 2 >F- £ " (, vH ) - 4v v + V 2’ 4 V >/ 2 ' I (4v v+ V 2 + 1 4 \>+ vj ) ' 


Thus, central differencing will do along the monotone decreasing profiles, and 
additional 0 ( |Av^ + y | ) amount of viscosity is required along the 
monotone increasing ones. 


3. FULLY-DISCRETE SCHEMES 

We consider two-level fully-discrete explicit or implicit schemes in the 
conservative form 


(3.1) v(t+At)=v(t)-X[h i.-h i.]. 

v v v v+ 72 v- V 2 

Here, Ax = i- (x . , - x , ) and At are the variable meshsize and time step 

v 2 v+1 v-1 

such that X y = » an< i h + is the consistent Lipschitz continuous 

numerical flux which depend on 2p+l neighboring gridvalues from both time 
levels, t and t + At. 

To study the TVD properties of these schemes, we forward difference (3.1) 


< 3 - 2 > 4 VV 2 <t+4t) 


4 v v + V2 U> ‘ Vl lh v + 3/2 ‘ Vl/2 1 + X v Ih v + V 2 ‘ V ‘/2 1 ’ 
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raultiply by s v+ l/^ + At) = sgn[Av^ + j^(t + At)] and sum by parts, 

obtaining 


(3.3) TV[v(t + A t) ] = 1 s v+ y(t + At) • Av tJ _ 1# _(t) + 


V W 


+ I V [ %+ l/ 2 (t + At > " s ,._ l/„ (t + At >l' f h ,.x l/„ - h .._ 1/J 


-1/ 2 V 


v+V 2 W 2 J 


The first summation on the right does not exceed TV[v(t)], and the requirement 
from the second one to be negative yields 


Leona 3.1: The fully-discrete scheme (3.1) is TVD, if we have 


(3.4a) 

h v+ V 2 — 

h v— l/ 2 

at maxima values 


+ At), 

(3.4b) 

h v+ V 2 — 

h v- I /2 

at minima values 

v v (t 

+ At). 


Lemma 3.1 is a manifestation of our previous assertion, namely, that the TVD 
properties of conservative schemes are determined solely by the behavior of 
their numerical fluxes at extrema values. Yet, unlike the semi-discrete case 
we had before, here there is the additional difficulty of tracing these 
unknown extrema values at the next time level, t + At. 

A similar situation occurs with the incremental representations of fully- 
discrete nonlinear schemes. Consider for example two-level explicit schemes 
in the incremental form 

v (t + At) = v (t) + C + . i. Av i,(t) - C 1 , Av i/(t), 
v v v+V2 v+ 72 v- 72 v-V2 ’ 


(3.5) 
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with C* = C* i,(v . (t)). To study the TVD properties of 

v+ 72 v+ V 2 v-p+1 v+p 

such schemes, we forward difference (3.5), multiply by s^ + y(t + At) and 
sum by parts, obtaining 

(3.6) TV [v(t + At)] = l s v+1/2 <t + At) • Av v+i/2 (t) + 

- I s v+ l/ 2 (t + At)s v + l/ 2 (t)[x v (t + At) C v+V 2 + X v+l (t + At)C v+V2 ] * ,AV v+V 2 (t)| * 

Since the first summation on the right does not exceed TV[v,(t)], we arrive at 
the following sufficient TVD condition. 

Lenma 3 . 2 : The explicit scheme (3.5) is TVD, if we have 

(3.7a) x v (t + at)C+ +I/2 + x v+1 (t + 4 C)c; + 1 /2 > 0, when V 2 (t + 40 " V 1 ^ 0, 

(3.7b) x v (t + ‘Oc; +1/2+ x v + ,(t + *t)c; +1/2 < 2. when « v+ + it) . ^(t). 

Applying the last result to the incremental coefficients, compare (2.28), 


(3.8) 


± 1 A f y+ V 2 

C V+ V 2 "2 ^ + V 2 ± Xa v+V 2 ' ) ’ a v+V 2 _ Av v+ 


we find that for equally spaced explicit schemes given in the viscosity form 



- 21 - 


(3.9) v v (t + it) - v v (t) - | [f(v v+1 ) - f( Vl )1 + 

+ 7 [ V V 2 1V v+ ‘/2 ' V V 2 ‘V W ' X v $ X ' 

the following TVD characterization holds. 

Lemma 3.3: The explicit scheme (3.9) is TVD, if its viscosity 

coefficient satisfies 


(3.10a) 


bc v (t 


+ At) +x v + 1 (t + At)].Q v+1/ 2 > [ Xy (t + At) - 


x v+l (t 


+ At) ] »Xa 


v+ V? 


and the following CFL-like condition is fulfilled 


(3.10b) X|a v+ i /2 |] < 1. 

Thus, we conclude that the TVD property of either scheme, (3.5) or (3.9), is 
determined by the behavior of their incremental and viscosity coefficients at 
extrema values, but as before, the difficulty lies in obtaining a priori 
knowledge about these values at time level t + At. 

The inequality 0 < x <2 suggests one way of avoiding this 
difficulty, namely, to replace (3.7) by the simpler positivity condition 

< 3 ' Ua > < + l/ 2 >°. C v + V 2 i°' 


together with a CFL restriction 



Yet, the simplicity of this sufficient TVD condition, which is originally due 
to Harten [3], [5], [9], [15], is obtained at the expense of its global 

dependence on the special incremental form being used. 

Another attractive approach to circumvent the difficulty of tracing the 
next-time level extrema values is to view the scheme (3.1) just as a first 
predictor step. Then, the resulting spatial variation at time-level t + At 
can be made the basis for an augmenting corrector step which will preserve the 
monotonicity of the predictor step and which will comply with (3.4). Such 
argument was used in connection with the FCT algorithm [1] and the ACM method 
[2]. In the following example, borrowed from [10, Corollary 4.9], we work out 
another corrective-type recipe of this kind, which highlights the essential 
features distinguishing the fully-discrete explicit case from the semi- 
discrete one. 

Example 3.4 : Consider an explicit first-order accurate TVD scheme of the 

form 

(3.12a) v*<t + 4t) - v v (t) - Mh u+1/2 - h v _ ,^1 - 

* V 0 - \ [f< Vn> • ‘'Vi 11 *7 l VV 2 4V v+V 2 " V 1/2 V 1/2 1 ' 

The star indicates the predicted values at time level t + At. In order to 
convert this scheme into a second-value accurate TVD one, these values are 
corrected to second order accuracy, by augmenting an anti-diffusive corrector 
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step of the form 


(3.12b) v v (t + At) = v y (t + At) - [^ v+1 / 2 " K v-1/ 2 ] * 


The numerical flux correction, Ti 1 , is chosen to be, compare Example 2.6, 

v+ 72 

«- 12 => Vv 2 -7 [S v + ?v + i - V#v + 1 - s.11- v*v + 


where should satisfy the two properties 


(i) sgn(g ) = s i, = s i, , at nonextrema values v (t + At) 
° °v v- V2 v+ 72 v 


(ii) g v = 0, 


at extrema values v^(t + At). 


In addition, we require that the predicted monotonicity should be preserved, 


i.e. , 


%+!,$<■ ♦“> - V‘/2 U + “)’ 


so that by the usual summation by parts we have 


TV[v(t + it)] - I s v+1/2 (t + it)i» v+1/2 (t + it) + 


* l U* +1/2 (t + it) - s v . 1/2 Ct + it)] . [h v+1/2 - h v . 1/2 1 


-24- 


Since the first summation on the right equals TV[v (t + At)] < TV[v(t)], 

while the second is nonnegative, consult Example 2.6, the scheme (3.12) is 
TVD. Next, its second order accuracy is achieved if 


\+ 1/ 2 *I [ W 2 " V V 2 ] ' Sv v+V 2 (t + 40 + 0[4v *+l/2 (t + 4t,]2 ‘ 


To satisfy this (away from extrema values), and the first two properties 
listed above, we choose 


(3 . 12d) g v -f V (Q v-V 2 - ‘V^’I'Vh! 1 


where the Lipschitz continuous form, B [ • , • ] , is yet to be determined so that 
the third property of raonotonicity preserving will be satisfied. To this end, 


we note that 




and therefore, since h i. and h i . must agree jlh sxgn, 

v+V2 v-V2 

l 4 v v+ V 2 (t +4t) ' 4 Cv 2 (t +4t)| ' l ’V3/2 - 2 V 1/2 + V 1/2 1 - iS vl + l*v+ll- 

Hence, the sum on the right does not exceed |Av^ + i^(t + A t ) | and 
consequently the monotonicity preserving property holds, provided B[»,«] 


must agree in sign, 


is chosen as the bilinear limiter form 


(3 . 12e) 


B[w 1 ,w 2 ] =Min[|w 1 |, |w 2 | ] . 


The last example demonstrates the typical situation with explicit 
schemes, where the TVD property necessitates one kind or another of a Minmod 
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limiter in order to prevent new extrema values other than those which 
propagate from time level t. 

An implicit version of the above corrective procedure is given in 


Example 3.5: Consider an implicit first order accurate TVD scheme of the 


form 

(3.13a) v*(t + At) + ^ [f(v* +1 ) " + 

+ i [Q i.Av " Q i,Av I,] = v (t). 
2 v+ 72 v+ 72 v- 72 v- 72 v 


We augment it with an antidif fusive fully-implicit corrector step of the form 


(3. 13b ) v v (t + At) + [\ + l/ 2 ~ \-l/} = + At >’ 

where 


(3.13c) 


V l/ 2 - 2 C % + «v + l ' V V 2 (t + 4t) ' l*v+l ’ 1 1 


Then (3.13b) serves as a second-order accurate solvable correction, if we set 


(3.13d) 



\<t + 4t).Minl(q vll4 + + »t)).|iv u±1/2 (t + 


The resulting scheme (3.13) is TVD under the original (possibly unlimited) CFL 
condition. Indeed, we have 
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TV [ v(t + At)] = l s v+1/ a + At)Av* +1/2 (t + At) + 

+ I [% + l/ 2 Ct + At) - s v _ 1/2 (t + At) ] • [h v+ 1/2 - h v _ 1/2 ] ; 

* 

the first summation on the right does not exceed TV[v (t + At)] _< TV[v(t)], 

<v 

while the second vanishes since h i, do at extrema values where 

v± V2 

s v+l/ 2 (t + At) * S v-V 2 (t + At) * 

Other recipes for constructing implicit TVD schemes which are second- 
order accurate away from extrema values, are suggested by the following 
analogue of Lemma 3.3. 

Lemma 3.6: The implicit scheme given in the viscosity form 

(3.1.4) v v (t + At) + ^ [f(v v + 1 (t + At)) - f(v v-1 (t + At))] + 

' 1 lq v + l/f v v. l/ 2 (t + 4t> ' V >/2 V l/ 2 (t + 4t>1 ‘ V v (t> 

is TVD, if we have 

% 

(3.i5) [X v (t + At) + x v+1 (t + At)]Q v+1/2 > 

[X v (t +At) -X v + 1 (t + At)]-Xa v+l/2 (t +At). 

We omit the proof and turn to our final 
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Example 3.7 : The viscosity of the second-order accurate implicit Lax- 

Wendroff scheme is modified at extrema values, by setting 


(3.16) 


1/2 


2 2 


= - X a 


+ '/a + 2 ( *v 


+ W 


[x la 


v+ 


vj + i 2 * 2 


v+ !/ 2 J 


J, 


where the quantities on the right are evaluated at time level t + At. The 
resulting scheme (3.14), (3.16) can be easily checked to satisfy (3.15) and 
hence is TVD. However, the linearized implicit LW scheme is unconditionally 
unstable — the amplification factors of its nonconstant modes, they all lie 
outside the unit disc [19]. Consequently, the TVD property of (3.16) is 
achieved by switching to upwind differencing at the extrema of these unstable 
oscillatory modes, at the expense of lowering the effective overall accuracy. 


Remarks : We note that the TVD characterization in Lemma 3.1 does not 
assume the CFL condition; it enters, indirectly, through the requirement of 
controlling extrema values at the next time-level. Substitution of the 
canonical incremental decomposition (2.14) into (3.7) reveals that the same is 
true with respect to the TVD conditions in Lemma 3.2 and 3.3, where the CFL 
limitation is implicitly contained already in (3.7a) and (3.10a). 
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